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Abstract 

Statistical properties of parametric motion in ensembles of Hermitian 
banded random matrices are studied. We analyze the distribution of level 
velocities and level curvatures as well as their correlation functions in the 
crossover regime between three universality classes. It is shown that the sta- 
tistical properties of level dynamics are in general non universal and strongly 
depend on the way in which the parametric dynamics is introduced. 



I. INTRODUCTION 



A link between random matrix theory (RMT) M and the statistical properties of spectra 
of quantum systems is well established. Depending on the symmetry of a classically chaotic 
quantum system, its spectral fluctuations are described by Gaussian orthogonal (GOE), 
Gaussian unitary (GUE) or Gaussian symplectic ensemble (GSE) 0,0]. 

Quite often the physical systems depend on some external parameter, say A, therefore it 
is interesting to study the level dynamics i.e. the motion of eigenvalues Ei(X) as a function 
of A. Among the first parametric properties studied were the investigations of the avoided 
crossings gaps JlJj]|§, the parametric number variance or the curvature of the levels (i.e. 
the second derivatives of their energies with respect to the parameter) |]^,|9|, [T^Jl~T| , p!2|Jl3|| . It 



has been claimed that the statistical properties of level dynamics are universal for disor- 
dered or strongly chaotic systems |8|,|l4}| provided the change of A does not modify global 
symmetries properties. To reveal the universality one has both to unfold the energy levels 
and appropriately rescale the parameter, A P, |IT|,[14| . Other statistical measures of paramet- 
ric dynamics such as the level slopes (velocities) distribution (Gaussian shaped for random 
systems p| JTT| , p^| ), the velocity- velocity correlation function []l4||15l , |16|,|17| in the bound spec- 



trum or parametric conductance fluctuations |T8[ and fluctuations in the Wigner time delay 
|H| for scattering systems have also been discussed. 

A word of caution is, however, necessary at this point. Even for the nearest neighbour 
spacing distribution, widely considered to be universal, exceptions from the RMT prediction 
may be quite significant for real physical systems f20 |. Much more pronounced and com- 
mon are the deviations from the RMT predictions for the parametric motion of levels. In 
particular, as shown by Takami and Hasegawa ||10|| , the curvature distribution shows nonuni- 
versal behaviour for small curvatures even for the mixing system (a Bunimovich stadium). 
Similarly non Gaussian slopes distribution as well as strongly non Cauchy-like curvature 
distribution was observed for magnetized hydrogen atom [fLlJ| . The origin of these devia- 
tions has been linked to partial wavefunction localization on unstable periodic orbits. Thus 
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the non-generic features of parametric statistics may provide most interesting information 
about the physics of a given physical system. 

Most of these studies considered pure symmetry cases, i.e., systems pertaining to a 
given, e.g., GOE universality class. This is often not the case in a realistic situation. In 
particular the presence of the magnetic field, or the Aharonov-Bohm flux the time-reversal 
invariance symmetry (TRI) becomes broken - such a situation corresponds to a crossover 
between GOE and GUE for random system. In this context the velocity-velocity correlation 
function has been studied intensively |2l|j2^JT7[1 as well as the velocity distribution |23] or the 



curvature distribution p4[^5|. The authors considered mostly the situation when the increase 
of the external parameter, A (e.g. the magnetic field) destroys the time-reversal invariance 
although, importantly, it has been noticed p3[ that the parametric velocity distribution may 



strongly depend on the nature of the perturbation. Relatively less frequent were studies 
of the parametric dynamics in the transition region between completely delocalized and 
localized spectrum (see, however, e.g., the treatment of the velocity distribution for broken- 
TRI case in [EBl 



In order to model the spectra of quantum systems in a crossover regime (a weak lo- 
calization or a partially broken symmetry) one may utilize random matrix ensembles that 
interpolate between the canonical ensembles. For example, real symmetric band random 
matrices are capable to model transition between the localized and the delocalized regime. 
Statistical properties of their eigenvalues and eigenvectors depend on a single scaling pa- 
rameter x = b 2 /N where iV denotes the matrix size and b the band width. Allowing 
the matrices to be Hermitian and changing the relative weight of the imaginary component 
a one can model the effect of the time-reversal symmetry breaking and the transition from 
an orthogonal to an unitary universality class. The corresponding scaling parameter y is 



proportional to Na [^] for small perturbations. An ensemble of Hermitian band random 
matrices (HBRM) can be therefore completely characterized by two scaling parameters (x, y) 
]. The Poissonian, strongly localized spectrum is obtained in the limit x « 1, while in 



the opposite delocalized limit x » 1 the model reduces to GOE for y — and GUE for 



y » i. 

This work is intended as a systematic study of the parametric dynamics and the cor- 
responding statistical measures for the transition region and in the localized regime. Our 
work differs from most of the analysis mentioned above in a way a parametric dependence 
is introduced. We assume that the changes of A leave the global properties of the system 
unaffected. In the random matrix approach this is equivalent to the assumption that the 
statistical properties of the ensemble of matrices do not depend on the value of the param- 
eter, A, determining the parametric dynamics. For the physical system applications this is 
equivalent to saying that the symmetry properties of the system considered, but also the 
character of the underlying classical dynamics (say, the fraction of the phase space volume 
which is chaotic) are invariant with respect to A. Or, from a practical point of view, the 
the dynamics changes only weakly with A, in the interval of A values considered in each 
case. The level dynamics is, in a sense, "perpendicular" to crossovers between canonical 
ensembles, as schematically shown in Fig. [I]. Such a physical situation may correspond to a 
variation of the disorder parameter in a mesoscopic system, for which all other parameters 
are kept constant. A special attention is drawn to the localized case, characterized by small 
values of x where some analytic predictions obtained using supersymmetric calculus exist 
2|g7|. 

The interest in such a study is twofold. Firstly, it is interesting to see how the parametric 
properties of the system follow the transition between different pure universality classes. 
Secondly, the results obtained in the HBRM model may serve in future as a reference for a 
comparison with statistics obtained in real physical systems. It is then of utmost importance 
to know what one may expect from the purely random model. This may enable to isolate 
the non generic, i.e. characteristic for a given system, properties. 

The paper is organized as follows. In the next section we describe the model and the 
parametric dynamics. The distribution of level velocities P(v) is analyzed in Section III. The 
next Section is devoted to velocity-velocity correlation function C V (X). Level curvatures are 
discussed is Sec. V, whereas their correlations are investigated in Section VI. Next Sections 
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consider higher order statistical measures. Finally we discuss the consequences of the results 
obtained for various statistical measures in the concluding section. 



II. PARAMETRIC DYNAMICS FOR HERMITIAN BAND RANDOM MATRICES 

Hermitian band random matrices are defined by 

By = tfS + i&Mb - \i - j\) ij = l,...N, (2.1) 

where 0(.) denotes the unit step function vanishing at the origin. Independent random 
variables and £/„■ are distributed according to Gaussian distributions with zero mean 
the root mean squares equal to af and afj, respectively. The parameter a measures the 
relative size of the imaginary part of the off diagonal matrix elements a = (ajj/a^) 2 , i ^ 
j; (notation has been simplified with respect to Ref. [30[|). A normalization condition, 
Tr (H 2 ) — N + 1, keeps all the eigenvalues in a constrained energy range. It also allows us 
to express the variances of real and imaginary parts of matrix elements in terms of matrix 
size N, integer band width b and real parameter a 



^= 2 N+{ JWK-W-D ll + Slih <2 ' 2) 



{(7ij) = 2N + (a + 1)(2N -b)(b- 1) (1 " 5ij) - (2 ' 3) 



For a diagonal random matrix (b = 1) the density of eigenvalues is Gaussian and the 
level spacings are distributed according to the Poisson distribution, independently of the 
parameter a. In the opposite limiting case of the full matrix (b = N) variations of the 
parameter a correspond to the process of the time reversal symmetry breaking in a dynamical 
system and control the transition between orthogonal (a = 0) and unitary (a = 1) ensembles. 

Statistical properties of spectrum and eigenvectors of real symmetric band matrices de- 
pend only on a scaling parameter x = b 2 /N. This scaling law, observed first by numerical 
computation of localization length [J2^J , was reported to describe also the distribution of 
eigenvalues ^TJ and eigenvectors |32| , and subsequently explained theoretically |33| . 



The same scaling law holds also for Hermitian matrices p3| , |34| . Moreover, effects of the 
time-reversal symmetry breaking are controlled by another scaling parameter y = 2Na/(l — 
a) [P0[], stemming from the universal properties of orthogonal-unitary transition founded by 



Pandey and Mehta ||29|| . The structure of eigenf unctions of HBRM and the distribution of 
inverse participation ratio has also been studied recently |5S| . 
Let us now consider the parametric random matrix 

H(X) =#i cos A + # 2 sin A. (2.4) 

Both matrices H\ and H2 are taken from the same ensemble of HBRM. Hence the spectral 
properties of H are stationary and do not depend on A. Moreover, during the transition, 
the motion of eigenvalues is restricted to a bounded energy interval for arbitrary A. This 
model of parametric dynamics was already used for investigation of level curvatures and 



velocity correlation functions |fl6| , p!7| . The dynamics of eigenvalues as a function of A may 
be treated as the dynamics of interacting particles (eigenvalues) with A playing the rule of 
the fictitious time 0. This allows to interpret the slope of the levels as the velocity of the 
particles and their curvature as the corresponding accelerations. 

Parametric dynamics defined above can be studied numerically is a straightforward way. 
For several values of ensemble parameters (N, b,a) we have generated random matrices 



according to Eqs. (|2.1| - |2.4j ). Diagonalizations of resulting matrices for several values of A 
have allowed us then to find level velocities and curvatures by a finite difference method. A 
special care has been paid to obtain reliable values of velocities and curvatures, especially for 
very small and very large values, by varying the size of the step in A [[3(J . Before computing 
the derivatives of eigenvalues with respect to A the standard unfolding technique was applied 
to set the mean level spacing A to unity. We have considered matrices of size N varying 
between 50 to 500, velocities and curvatures have been computed at about 200 different 
values of A, and the typical number of the independent realizations of dynamics, Eq. ( f2.4| ). 
in each case studied has varied with matrix size to ensure at least 200 000 data in each 
statistics. In other words, we have simultaneously performed the averaging over the energy 



(data from different energy levels of a given matrix H) and the averaging over the disorder 
parameter (several realizations of the dynamics for the same values of N, b and a). 

To check reliability of the numerical procedure we have set b = N and we have reproduced 
known results concerning the distribution of velocities and curvatures as well as the velocity 
correlation function for GOE (a = 0) and GUE (a = 1). Moreover, we have verified 
that both scaling parameters, x and y, correctly describe the parametric dynamics. The 
statistical properties of all quantities studied have been found to be independent of the 
matrix dimension N (for sufficiently large N) provided the parameters x and y have been 
kept constant. 

In the following sections we describe novel results obtained for different statistics, com- 
mencing with the distribution of first derivatives, i.e. velocities. To avoid any misunder- 
standing let us repeat again that all the data presented are obtained for "perpendicular" 
transitions [both Hi and H 2 in Eq. (|2.4j) belong to the same random matrix ensemble] as 
exemplified by double-sided arrows in Fig. |I[ Thus for all values of A the scaling parameters 
x and y have the same values. We shall not consider here the case when the parameter 
change modifies the global symmetry properties - a situation exemplified by broken line 
arrow in Fig. [I]. 



III. DISTRIBUTION OF LEVEL VELOCITIES 

For level dynamics within GOE or GUE the distribution of level velocities, P(v), is 
Gaussian |^, p] , ^| . This fact is easy to explain using the first order perturbation theory. 
For A = the derivative dEi/dX is equal to the diagonal element of matrix H 2 expanded in 
the eigenbasis of Hi. Since both matrices are drawn independently from the same ensemble, 
the matrix elements are Gaussian random numbers leading to Gaussian velocity distribution. 

On the other hand, in strongly localized limit an analytical formula for P(v) given by Fy- 



odorov p6[ for systems with a broken TRI strongly differs from a Gaussian. A non-Gaussian 



character of the velocity distribution for the GOE — > GUE transition, corresponding to the 
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TRI symmetry breaking, has been discussed in p3| , |2^j25[ . 

We have analyzed the transition between localized and delocalized spectra both for ran- 
dom systems with a broken TRI (i.e., the ensembles interpolating between Poisson and 
GUE) and for ensembles interpolating between Poisson and GOE. The former allows us to 
test the analytical prediction of Fyodorov [EH] . 



The theoretical prediction, as presented in p6| , has no free parameters, both the shape 



of the distribution and its scale (determined by the velocity variance) are determined by 
the theory. Surprisingly the direct comparison of that distribution with the numerical data 
obtained has been highly unsatisfactory. The agreement is recovered, see Fig. |2|, when both 
the theoretical distribution and the numerically obtained data are rescaled with respect to 



the velocity variance, a v = y {v 2 } (note that the mean velocity vanishes by the construction 
of the ensemble) . Thus the apparent disagreement originally observed is due to the difference 
between the theoretical and numerically obtained values of the velocity variance (the ratio 
of the numerical value to the theoretical prediction being about 13). We do not have a 
clear explanation of this disagreement. It may be due to the fact that the band width in 
our HBRM ensemble is sharply defined [compare Eq. ( |2.1| )] while Fyodorov [26] assumed a 



smooth decrease of the random matrix elements variance with increasing distance from the 
diagonal, \i — j\. 

The theoretical prediction |26| , represented by a smooth line in Fig. 0, takes the form 

7T 7TWCOth(7Tw/v / 6) - V 7 ^ . < 

P ^ = ft r-T27 T7W\ ' ( 3 - 1 ) 

o smh (7TW/V6) 

where the rescaled velocity w = v/a v . The similar quality agreement is obtained for different 
values of the scaling parameter up to x of the order of unity corresponding to the transition 
to a delocalized case. Then the numerical data start to show Gaussian (typical for GUE) 
large velocity tail instead of the exponential tail corresponding to fully localized situation, 
as exemplified by the lack of large velocities in the numerical data presented in Fig. EJ as a 
thin line histogram. 

Although the theoretical prediction is obtained for the case of a fully broken TRI our 
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numerical data indicate that it works extremely well also for preserved TRI (real symmetric 
matrices) provided that again the velocity variance is appropriately adjusted. The exemplary 
data are presented in Fig. [3] for two cases corresponding to strong localization and a transition 
to the delocalized regime. Here the numerically obtained variance is twice larger than the 
theoretical value calculated in the same ways as for the broken-TRI ensemble. It seems, 
therefore, that the same distribution, Eq. ( jOl) , describes the velocity distribution for both 
TRI case and the no-TRI situation. The difference between the two ensembles (the former 
interpolating between Poisson and GOE, the latter between Poisson and GUE) appears 
in the numerical value of the velocity variance only. It is clear the variance is a unique 
parameter that determines the appropriate velocity scale, similarly as for GOE and GUE 
pure ensembles ||. 



IV. DISTRIBUTION OF LEVEL CURVATURES 

Let us consider now the distribution of curvatures, K = d 2 E/d\ 2 . As shown by Gaspard 
and coworkers || the tail of the distribution decays algebraically as K^ 2 ^ 13 . This universality 
has been verified for different systems P JT0|jn| ]. At the same time, the small curvature 
behaviour has been found to be non-generic even for strongly chaotic systems []TTJ],[TT|] and 
reflecting the system-dependent wavefunction localization properties (scarring by periodic 
orbits). 

On the other hand, the scaled curvature 



A 



for pure random ensembles obeys the generalized Cauchy distribution [11,12.13 



P{k) = N P - * , (4.2) 

(1 + K l ) 2 

(with p = 1, 2, 4 for GOE, GUE, GSE, respectively, Np denotes the normalization constant). 

Here we demonstrate that if one allows the parameter (3 to acquire real values, (3 G (0, 2], 
the same distribution may be used in a general case of the intermediate ensemble interpolat- 
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ing between Poisson, GOE and GUE pure cases (provided we consider the "perpendicular" 
transition). The normalization constant is then equal 

and the rescaling, Eq. ( |4.1| ) holds almost everywhere. 

To test this conjecture we have generated several parametric dynamics "perpendicular" 
to crossovers between Poisson-GUE, Poisson-GOE and GOE-GUE using, as before, the 
formulation of Section II, Eq. (|2.1| ) and Eq. ( |2.4j) . The numerically obtained histograms of 
curvatures in the double logarithmic scale have been used to fit the algebraic decay of the 
tail of the distribution to the formula P(K) ~ . Then f3 has been found as (3 = fi — 2 
[compare Eq. (|4.2j )1. The same value of (3 has been used, together with the numerically 
obtained velocity variance to rescale the curvatures according to Eq. (|4.1|) . The exemplary 
results of such a procedure together with the conjecture ( |4.2|) are presented in Fig. [| and 
Fig. |5| in double linear and double logarithmic scale, respectively. Observe the excellent 
agreement between the numerical results and the proposed distribution. 

While Eq. (|4.2|) seems to describe well, at least approximately, the numerical data for 
curvatures everywhere in between pure cases of GOE, GUE and Poisson limit, the scaling 
Q4.1D works best for the delocalized or weakly localized spectra. For the Poisson - GUE 
crossover, close to the Poisson limit, the scaling obtained using Eq. (|4.1| ) is incorrect. The 



agreement with the generalized Cauchy distribution (|4.2j ) is obtained only if the numerical 
data are rescaled additionally by a numerical factor of the order of unity (typically 1.5 - 
2, depending on (3). Putting it differently, the parameter (3 in the denominator of Eq. ( [4. 1| ) 
should be replaced by other function of (3 which goes to (3 when transition to delocalized TRI 
broken case (i.e. GUE) is fully accomplished. This indicates that the proposed distribution 
( |4.2| ) is most probably the approximate one only. Still we find it quite remarkable that this 
simple analytic expression, with the proper rescaling, works so well for the interpolating 
ensembles. 

Let us mention that the power of the algebraic tail behaviour may be analytically related 
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to the level repulsion parameter (3 by a simple consideration of 2 x 2 random matrices || 
yielding fi = (3 + 2. A comparison of (3 values obtained from the tail of the distribution 



with (3' values obtained from the independent fit of the Izrailev distribution Pp>(s) [53] is 
presented in Fig. ^| in the whole interval of the intermediate f3 values. The agreement is 
quite good (and of the similar quality as that obtained for the Fourier transform of the 
velocity- velocity correlation function) considering that both the spacing distribution P/3>{s) 
and the proposed curvature distribution P(k) are most probably the good approximations 
to the true distributions only. 

It is worth noting that the distribution (|4.2| ) works well for the ensemble interpolating 
between GOE and GUE for the perpendicular action of the parameter A. On the other 
hand, if A is responsible for the TRI symmetry breaking, it has been shown that the tails 



of the curvature distribution are exponential [24,25 and not algebraic, as observed in this 



work. Parametric statistics are, therefore, sensitive to the way the parameter acts. Another 
example of this sensitivity are available from the earlier studies of periodic band random 
matrices [37| and 3-D Anderson model [39|, where the curvature distribution close to a log- 



norm distribution has been observed in the localized case while yet another distribution has 
been proposed in transition regime P0fl . 

V. VELOCITY VELOCITY CORRELATION FUNCTION 



In a series of papers, Altshuler and co-workers |L4| have discussed the universality of 
parametric statistical properties for disordered samples as well as for Gaussian ensembles. 
To reveal the universality both the eigenvalues (to unit mean spacing) and the parameter A 



(as X = a v X) have to be rescaled |8| JTT|JT4]1 . We have observed the power of such a rescaling 
already in the previous Sections. 

Consider next the velocity-velocity correlation function, a frequent subject of recent 
investigations [0,0,0,0,0,0,0, 



11 



1 2n 

C«(A) := / ^(A')^(A' + A)dY), (5.1) 



2ttA 2 



o 



where ( ) denotes ensemble averaging and A stands for the mean level spacing. By definition 
C v (0) = al thus the appropriately rescaled correlation functions take the form C V (X) = 
C V (X) I '<j%. Moreover, several models of time reversal symmetry breaking due to Aharonov - 
Bohm flux lead to correlation function practically indistinguishable from this characteristic 
for GUE. 

It was shown JT^ , pT| , |Tol] that for all three universality classes the rescaled correlation 



function 

C V (X) ~ A P X~ 2 , X -> oo (5.2) 

with the proportionality coefficient Ap dependent on the ensemble (we denote by ft the level 
repulsion parameter, ft = 1,2,4 for GOE, GUE and GSE, respectively). 



Explicit expressions have been obtained [14 for a closely related [but distinct from C V (X)] 



autocorrelation functions at fixed energy. A global approximation for C V (X) has been pro- 
posed |]IE|| . For the case of a classically chaotic system subject to a Aharonov-Bohm flux 
Berry and Keating |22| obtained a semiclassical approximation for C V (X) having the form 
of an everywhere analytic function of A. Yet it was demonstrated |T7[ that C V (X) is not 
analytic and suffers a logarithmic singularity at A = 0. 

Analytic properties of correlation functions are conveniently studied using the periodicity 
in A. In the Fourier domain 

oo 

C V {X) = ^2c n cos(nA). (5.3) 

ra=0 

Mean squared velocity, determining the scale, is given by the sum of all coefficients o 2 = 

oo 

C(0) = J2 c n- Expanding the dependence of a given eigenvalue on A in the Fourier series 

n=0 

Ei — S^L-oo a i emA ; where a n = a*_ n , on account of Eq. (|5.1| ) it is easy to see that c n = 
n 2 (\a n \ 2 )/A 2 . 
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Asymptotic behavior of C V (X) ~ — A~ 2 corresponds to a linear raise of Fourier coefficients 
c n for small n. On the other hand, for large n it was shown [IT] that (|a n | 2 ) ~ rT A ~^ and 
consequently, 

c n ~ n~ 2 -P. (5.4) 

This result was obtained extending the parameter A into the complex plane and analyzing 
the distribution of branch points and anticrossings (||| . 

Thus, despite nonanalytic character of the correlation function both C V (X) and its Fourier 
transform have simple asymptotics given by Eq. (|5.2p and Eq. (|5.4j) , respectively. It is 
interesting to see whether a similar behaviour may be found for the interpolating ensembles. 
To this end we have studied the asymptotics of C v for all three possible transitions, i.e., 
ensembles interpolating between Poison and GOE, Poisson and GUE, as well as GOE and 
GUE. In all cases the parameter A acted perpendicularly to a given transition (compare 
Fig. 0). 

We have observed the same X~ 2 large X behaviour, Eq. ( |5.2| ) independently of the 
ensemble studied. As an example Fig. |7J shows the rescaled velocity-velocity correlation 
functions C V (X) corresponding to 5 different cases along the Poisson-GUE crossover. 

Algebraic decay of the corresponding Fourier transforms is visualized in Fig. Observe 
the continuous change of the slope, growing from —4 for GUE till —2 for the Poisson limit. 
This corresponds to the continuous change of the repulsion parameter f3 between 2 and 
in the level spacing distribution. Independently one may fit the nearest neighbour spacing 
distribution obtained numerically to the Izrailev distribution [38[|, P ( g/(s), which provides an 



excellent approximation for the nearest neighbour spacing distribution for the interpolating 
ensembles. P/3'(s) ~ s@ for small s. We have verified that (3 values obtained by fitting 
the straight line to the tail of log(c n ) equal (3' values obtained from the fits of the spacing 
distribution within 5%. Therefore, we conclude that the validity of Eq. (|5.2| ) and Eq. ( p.4| ) 
extends to the intermediate ensembles and fractional values of the repulsion parameter (3. 
We stress again that this result is restricted to the "perpendicular" transitions only. 
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If parameter A is responsible for the transition between the ensembles (e.g., the magnetic 
flux in the Aharonov - Bohm effect) the velocity-velocity correlation function obeys the 
c n ~ n~ 4 algebraic decay, independently of the degree of localization [ IT] . Thus, similarly 
to the velocity distribution itself also the velocity- velocity correlation function C V (X) is 
sensitive to the nature of perturbation generating the parametric dynamics. 



In a full analogy with the velocity-velocity correlation function (|5.1| ) we define the cur- 
vature correlation function 

1 2w 

Ck(X) := ^2 (/ Ki(X')Ki(X' + X)dX'). (5.5) 
o 

However, this function does not provide us with any new information. This fact is easy to 
understand studying the Fourier expansion Ck(X) = J2^Lo k n e mX . As for velocity correlation 
function one uses mean Fourier coefficients of individual energy levels and obtains relation 
k n = n 4 (|a n | 2 )/A 2 . A comparison with the velocity correlation function Fourier coefficients 
yields immediately 

C k {\) = g^C v (X), (5.6) 

which easily yields the properties of Ck{\) from known properties (e.g. the asymptotic 
behavior) of C„(A). Eq. ( |5.6j ) holds for an arbitrary matrix ensemble. For completeness 
we present the numerically obtained Cfc(A), rescaled with respect to Cfe(O), for GOE and 
GUE ensembles in Fig. |9|. Notice a much faster decay of the correlation between curvatures 
as compared with the velocity correlation function. Asymptotically, using C v ~ A -2 and 
Eq. ( |5.6|) we get oc A -4 for large rescaled parameter A, in full agreement with Fig. ||. The 



curvature correlation function, in view of Eq. ( |5.6|) , may be used, together with the velocity 



correlation function, for numerical tests of the accuracy of curvature evaluation (which may 
be quite tricky using the finite difference method since small and large curvatures may 
require different step in the parameter). 
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VI. HIGHER DERIVATIVES OF ENERGY WITH RESPECT TO THE 

PARAMETER 



Algebraic decay of the Fourier transform of velocity correlation function on one hand 
provides an information about singularity of some higher derivative of C V (X) at A = 0, and 
on the other hand, indicates a possibility of a divergence of a distribution variance of the 
some higher derivatives of the energy levels with respect to the parameter E7| . In particular, 
for orthogonal ensemble (/3 = 1), the variance of curvature distribution (K 2 ) does not exist 
and in order to characterize the mean curvature one uses the mean absolute value (\K\) 
instead |37j . Moreover, the second moment of the distribution of third derivatives of energy 
levels L := d 3 E/d\ 3 was predicted |17| to diverge for (3=1 and (3 = 2. 

To test this prediction we have studied the distribution of these third derivatives. The 
most difficult part here is to find how to call them - using the level motion picture where the 
curvature of the level is identified with the acceleration of the fictitious particle, the third 
derivative of the energy will correspond to the derivative of the acceleration. In spirit of this 
mechanical analogy we refer to the third derivative as a jerk (editorial - referee's - reader's 
help here, whichever comes first, would be of great value for the authors). 

We have restricted the numerical study of the distribution of jerks to canonical orthogonal 
and unitary ensembles (GOE and GUE). The obtained numerical results are displayed in 



Fig. As expected the distributions of jerks are characterized by the algebraic tails; the 
numerically obtained power law decay yields P(L) ~ L^ +3 ' 2 ' . That confirms the divergence 
of the variances both for GOE and GUE. 

It is interesting that the distribution of jerks may be quite nicely approximated by a very 
simple ansatz 

P{L) = Np gas , (6.1) 

(1 + B^))WJ 

where Ma is a normalization constant while Bp and A(j3) are free parameters. In Eq. ( |6.1|) 
the jerks are conveniently rescaled taking the unfolded spectrum (i.e. with the mean level 
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spacing equal to unity) with derivatives calculated with respect to the rescaled parameter, 

x = ^F)\. 

We have fitted the distribution of this form to the numerical data for both GOE and 
GUE. The results are represented in Fig. |H] as smooth curves and quite successfully represent 
the numerical data. The obtained values of the parameters are equal for GOE to A\ = 1.67 
and B\ = 9.08 while for GUE we obtain A 2 = 2.50 and B 2 = 0.84. The obtained values 
of Ap are close to simple fractions, A± = 5/3 and A 2 = 5/2 - the corresponding curves are 
indistinguishable from best fits within the accuracy of our data. 



VII. CONCLUDING REMARKS 

We have analyzed various aspects of parametric dynamics in the space of Hermitian 
random matrices. Such a model may be applied to study transitions between Poissonian, 
orthogonal and unitary universality classes. We have analyzed numerically the situations 
when the parameter change does not modify the global properties of the ensemble studied, 
the case baptized as a "perpendicular" transition to contrast it with the "parallel" case 
when the parameter is responsible for the break up of the symmetry or other change of the 
properties of the ensemble studied. We have, however, compared our results with predictions 
of other works where often such a "parallel" parameter action was considered. 

We have paid a particular attention to the study of the transition between the Poissonian 
ensemble characterized by strongly localized wavefunctions and the delocalized Gaussian 
ensembles (GOE or GUE). In particular, the numerical tests of the analytic predictions for 
the distribution of level velocities in the case of broken time reversal invariance |26[ confirmed 



the predicted shape. We have observed a disagreement between the theory p6| and the 



numerical data, however, as far as the prediction for the velocity variance is considered. We 
have discussed the possible origin of this difference. We have shown that the distribution of 
the same functional form works for the derealization transition also for the real symmetric 
random matrices. This calls for the extension of the theory to such 
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We have studied in detail also the distribution of level curvatures. Here no analytic pre- 
diction is available. We have found that the numerically obtained distribution of curvatures 
is well approximated by the generalized Cauchy distribution (|4.2|) earlier shown to be exact 



TT| , [T2|JT3|1 for the canonical, GOE and GUE ensembles. The only required modification is to 
take the fractional value of the level repulsion parameter f3 in accordance with the spacing 
distribution. We have found also that the same rescaling ( |4.1|) holds everywhere except 
for the localized, no-TRI ensemble. Then the agreement with Eq. ( }4.2|) requires additional 
multiplication of all curvatures (rescaling) by a factor of the order of unity and dependent 
on f3. 

This form of the curvature distribution implies its algebraic tails of the form P{K) ~ 
j^-i-p Similarly we have found that the same level repulsion parameter f3 governs the 
tails of the Fourier transformed velocity correlation functions. Explicitly, the corresponding 
Fourier coefficients satisfy to a good precision c„ ~ rT 2 ^^ . 

Comparison with other works, where mostly the "parallel" transition have been studied 
PB| , |2^p5| , |5B]P7j , |55[ ] on various models indicates strong differences with the "parallel" tran- 
sitions. This difference has been first observed for the velocity correlation function in the 
case of a partially broken TRI in a fully delocalized case Here we have shown that the 
sensitivity of level dynamics to the way in which the parameter acts extends also to other 
parametric statistical measures as well as to other ensembles interpolating between "pure" 
cases of the Poisson, GOE and GUE ensembles. This has an important consequence — it 
shows that the universality of parametric dynamics is more limited that anticipated before 

0- 

Finally consider the consequence of the presented results for studies of realistic systems. 
Consider the semiclassical limit when the system is "large", with a high density of states and 
many highly excited levels. In the generic situation a small change of the parameter cannot 
induce significant changes in system symmetries and global properties - small changes of a 
parameter may be thus considered as "perpendicular" cases. This indicates that the "per- 
pendicular" transitions studied in this paper are typically generic. An important exception 
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is the change of the magnetic fields in systems with no additional symmetries where the field 
induces the breakup of the time reversal symmetry acting, therefore, in a "parallel" way. 

Thus the fact that the curvature distribution for the transition between the Poisson en- 
semble and GOE, realized via the banded matrix model studied here, is described by the 
generalized Cauchy law, Eq. (|4.2| ) has important consequences. It has been suggested (see 
pOfl and references cited there) that banded matrices may be used to simulate statistical 
properties of partially chaotic systems interpolating between the integrable case (with Pois- 
son level spacing statistics) and fully chaotic case (GOE). This has been partially based on 
the similarity of the level spacing distribution observed in both cases, well approximated 



by the Izrailev ansatz or for TRI systems by the Brody distribution. However, even 
for chaotic systems, as shown before, small curvature behaviour may be abundant due to 



isolated avoided crossings and scarring of wavefunctions |TT]. For the mixed phase space 
systems the avoided crossings are typically quite narrow and isolated - between them the 
levels can be adiabatically followed as a parameter is varied. Generically, small changes 
of a parameter are accompanied by small changes of eigenenergies which may be treated 
by a Taylor series expansion with a leading linear term. It shows that such systems will 
exhibit a great abundance of small curvatures, accompanied possibly by a singularity of the 



distribution at K = (see also for an additional discussion and numerical examples). 
Banded random matrices show a different curvature distribution than expected for a quan- 
tum system with a mixed phase space. Thus this ensemble is not adequate for simulating 
the statistical properties of partially integrable or weakly chaotic systems at least in cases 
when parametric dynamics is concerned. On the other hand, HBRM ensemble seems to 
be very useful for obtaining predictions for random systems which exhibit a transition to 
localization, well into the localization regime. 
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FIGURES 

FIG. 1. Scheme of the space of random matrices (and dynamical systems). Three circles repre- 
sent universality classes: Poissonian, orthogonal and unitary, and dashed lines represent crossovers 
between them. Solid arrows stand for "perpendicular" transitions, analyzed in this paper. A 
broken arrow exemplifies a "parallel" transition, not treated here. 

FIG. 2. Velocity distribution in a semilogarithmic scale for HBRM model interpolating between 
Poisson ensemble and GUE (a = 1, a fully broken TRI). Thick (thin) line histogram corresponds 
to numerical data obtained in the localized case x = 0.126 (the delocalized case, x = 1.408) from 
diagonalizations of 10000 matrices of rank N = 71. A dashed thick line represents the theoretical 
prediction, Eq. ( jSP| ) while a thin line represents the Gaussian distribution. 

FIG. 3. Same as Fig. 2 but for real symmetric matrices, a = 0: thick line histogram - a localized 
case, x = 0.056, while thin line histogram - a delocalized case, x = 5.63. 

FIG. 4. Exemplary curvature distribution for the ensemble interpolating between the localized 
Poisson case {(3 = 0) and GOE (/? = 1) obtained numerically for matrices of rank N = 71 and 
bandwidth 6 = 5 corresponding to x = 0.352 (histogram). Thick dashed line represents the 
fitted distribution, Eq. ( ^4.2| ) with (3 = 0.59. Thin solid and dotted lines represent the limiting 
distributions for the Poisson ensemble and the GOE, respectively. 



FIG. 5. Same as Fig. || but in a double logarithmic scale. 

FIG. 6. Parameter f3 obtained from the decay of the tail of the curvature distribution against 
the level repulsion parameter j3' obtained from the independent fit of Izrailev distribution [38]. 
Each dot represents one ensemble interpolating between the Poisson ensemble, GOE or GUE. 
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FIG. 7. Rescaled velocity correlation function C V {X) obtained for 5 transitions perpendicular 
to the GUE-Poisson crossover: N = 71, c = 2.0; b = 71(A) (GUE), b = 10(+), b = 7(a), b = 5(o) 
and b = 4(0). Ensemble averaging performed over 100 matrices; the lines are drawn to guide the 
eye. 

FIG. 8. Fourier transform of velocity correlation functions displayed in Fig. [7| in the log-log 
scale. Lines represent slopes characteristic for GUE (—4), GOE (—3), and Poisson (—2). The data 
for b = 10 are not plotted to improve legibility of the figure. 

FIG. 9. Curvature-curvature correlation function Ck(X) obtained for N = 50 (O), 60(o), 
70(Q), 80(A), and N = 90(<?/*) for GOE (open symbols) and GUE (full symbols). Universal 
rescaled velocity correlations are represented for comparison by thick dashed (GOE) and solid 
(GUE) lines. 

FIG. 10. Distribution of jerks (third derivatives of the eigenenergy with respect to the parame- 
ter) P(L) in the double logarithmic scale. Histograms correspond to numerically obtained data for 
GOE (thick line) and GUE (thin line) both obtained for random matrices of rank N = 71. Dotted 
and dashed lines represent the best fit of the proposed distribution, Eq. for GOE and GUE, 

respectively. 
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